DUDLEY KNOX LIBRARY 

W*f p^^'DR'Yifsrr */*M*gj 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




GENERATION AND ANALYSIS OF 
PARAMETERIZED TERRAIN FOR LAND COMBAT 

by 

Chung Young Kim 
March 1977 



Thesis Advisor: 



S. H. Parry 



Approved for public release; distribution unlimited. 



T178062 



UNCLASSIFIED 

SECURITY CLASSIFICATION OF this PAGE (When Dote Bntmrrnd) 



rm&fmwm 

hAVAl POSTGRADUATE SCHOOL 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


fTRERORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE 'ond Submit) 

Generation and Analysis of 
Parameterized Terrain for Land Combat 


5. TYPE OF REPORT A PERIOD COVEREO 

Master's Thesis; 

March 1977 


8 . PERFORMING ORG. REPORT NUMBER 


7. »uTh04(«; 

Chung Young Kim 


8 . CONTRACT OR GRANT NUMlERfij 


». performing organization name ano address 

Naval Postgraduate School 
Monterey, California 93940 


10. PROGRAM ELEMENT. PROJECT, TASK 
AREA A WORK UNIT NUMBERS 


11. CONTROLLING OFFICE NAME ANO AODRESS 

Naval Postgraduate School 
Monterey, California 93940 


12. REPORT DATE 

March 1977 


11. NUMBER OF PAGES 

59 


T A MONITORING AGENCY NAME A ADO RCSS(lf different from Controlling Office) 


IS. SECURITY CLASS, (ol MJ* report; 

Unclassified 


IS«. OECL ASSIFICATION/ DOWN GRADING 
SCHEDULE 


16. DISTRIBUTION STATEMENT (ol thlt Report; 

Approved for public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (of the mb m tract entered In Block 30, tt different from Reptrrl) 


is. supplementary NOTES 


19. KEY WORDS (Continue on rororoo mldo If nocommm y mnd Identity by block number) 

Parametric Terrain 


20. ABSTRACT (Continue on referee mldo If neceoemy end Identity by block number) 

Many current army combat simulation models require 
digitized terrain inputs for their execution, resulting 
in a costly and time-consuming operation. A methodology 
to generate parametric terrain as a function of the number 
of hill masses and the height, location, spread, and slope 
' of the hills is presented in this paper. Modified Bivariate 
Normal (MBVN) distribution function and a linear function 



DD 1473 EOITION OF I NOV 8* IS OBSOLETE 

(Page 1) S/N 0102-014-6601 | 



UNCLASSIFIED 

SECURITY CLASSIFICATION OF THIS RAGE (When Date Entered) 



1 



UNCLASSIFIED 

teCUHlTv CLASSIFICATION OF THIS P«0£|-<F>i«n D»t« Enf.r.J 



(20. ABSTRACT Continued) 

cutting the MBVN function are employed to achieve the terrain 
representation. Analysis of various terrain parameter values 
was conducted and the significant trends resulting from the 
study are described. Finally, recommendations for future 
enhancement of the model are given. 



DD Form 1473 
, 1 Jan 73 

S/N 0102-014-6601 



UNCLASSIFIED 

SECURITY CLASSIFICATION of this PAGErWh#n D»tm Enfrmd) 



Approved for public release; distribution unlimited 



Generation and Analysis of 
Parameterized Terrain for Land Combat 



by 

Chung Young Kim 

Lieutenant Colonel, Korean Army 
B.S., ROK Military Academy, 1963 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN OPERATIONS RESEARCH 



from the 

NAVAL POSTGRADUATE SCHOOL 
March 1977 



ABSTRACT 



Many current army combat simulation models require 
digitized terrain inputs for their execution, resulting 
in a costly and time-consuming operation. A methodology 
to generate parametric terrain as a function of the 
number of hill masses and the height, location, spread, 
and slope of the hills is presented in this paper. Modi- 
fied bivariate normal (MBVN) distribution function and a 
linear function cutting the MBVN function are employed to 
achieve the terrain representation. Analysis of various 
terrain parameter values was conducted and the significant 
trends resulting from the study are described. Finally, 
recommendations for future enhancement of the model are 
given. 
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I . INTRODUCTION 



Many current army combat simulation models require 
digitized terrain inputs for their execution. This process 
is very costly and requires a significant amount of computer 
storage. In addition, computation of line-of-sight and 
mobility require a great deal of computer running time. 

The current methodology utilized to represent terrain 
is also questionable from a statistical point of view. 

A Vector Research Incorporated (VRI) study on terrain line 
of sight, concluded that "... present and past army study 
results, based on the analysis of combat results on a very 
limited sample of terrains, may have been determined by 
the terrain selection process and not by the actual weapon 
system or force design differences." More specifically 
it stated the following: 

A. "There is extreme sensitivity in combat model 
results as the scenarios (terrain and movement assumption) 
are varied, even when variation is within a class of 
scenarios chosen for their a priori equivalence." 

B. "This sensitivity can be slightly reduced, but 

remains extreme (with probabilities of win estimable only 

within plus or minus 25%) even when battle results are 

„1 

used to redesign scenarios." 

^"Farrell, Robert L. and Richard J. Freedman, Investigation 
of the Variation of Combat Model Predictions with Terrain Line 
of Sight, Vector Research Incorporated, 1975. 
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These results imply that sufficient replications of each 
type of terrain should be run in order to reach a satisfactory 
statistical level of significance. For VRI's analysis at 
least fifty replications of each type of terrain were used. 

In view of the problems stated above, research was under- 
taken by Major Christopher Needels to develop a new method- 
ology for the representation of terrain. In particular, a 
modified bivariate normal distribution function was utilized 
to generate hill masses in a wide variety of configurations. 
Line of sight and movement routines were also developed to 
evaluate various intervisibility parameters over many terrain 
configurations for very low cost. 

This thesis represents an extension of the models 
developed by Needels to generalize the terrain generation 
routines. In particular, methodology was developed to 
generalize the hill mass generation from symmetric hill 
shapes to non-symmetric configurations. The methodology 
developed for this generalization is described in Chapter III. 

In order to evaluate the effects of the number of hills, 
hill height, and slope on various intervisibility parameters, 
an experiment was designed to measure these effects. A 
description of the design methodology and results of the 
analysis are given in Chapter IV. Finally, the conclusions 
drawn from the current analysis and recommendations for 
future research are presented in Chapter V. 

2 Needels, Christopher J. , Parameterization of Terrain in 
Army Combat Analysis , Master's Thesis, Naval Postgraduate 
School, Monterey, California, 1976. 
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II. REPRESENTATION OF TERRAIN 



A. HIGH RESOLUTION COMBAT MODEL REPRESENTATION 

The classic approach to representing terrain in land 
combat models is to utilize digitized terrain evaluation 
data from Waterway Experiment Station (WES) for the 
particular area of interest. 

The Dynamic Tactical Simulator (DYNTACS) model serves 
as a good example of the state of the art in terrain 
representation and use. DYNTACS is a two-sided dynamic 
model of battalion level combat representing the detailed 
interaction of elements in a combined arms environment. 

DYNTACS utilizes 100 meter grid squares for which 
macro-terrain elevations are determined from WES tapes . 

The simulation divides each grid square diagonally, thus 
providing a series of adjoining triangular terrain. The 
entire battlefield, therefore, is represented as a surface 
of equal size triangles which vary in slope depending on 
the elevation at their corners. 

Needels (pp. 9-10) provides a concise description of 
the DYNTACS methodology for line of sight calculation as 
follows . 

From the macro-terrain data, line of sight 
between any two opposing elements is computed. 

This is accomplished by first computing the 
angle between the horizontal and a straight 
line drawn between an observer and target 
(O-T Line) . The program then conducts a 
search of the terrain along the path of the 
O-T line to see if any macro-terrain is 
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higher than the 0-T line itself. This is 
accomplished by comparing the angle of the 
0-T line with the angle above horizontal 
of the Observer-Terrain line. If the latter 
angle is larger, there is no intervisibility. 
Over the duration of a battle with numerous 
elements this calculation may be made 
thousands of times. Consequently, not only 
is the time to prepare the tapes high, but 
also the time to compute lines of sight once 
the terrain is input to the model. 
Considerable effort by developers and users 
of this model has been expended to stream- 
line this subroutine. 



B. SYMMETRIC REPRESENTATION 

3 

The basic models developed by Needels (hereafter referred 
to as SIMTER) which provide the basic foundation for the 
methodology developed in this thesis are described in this 
section. The basic motivation for SIMTER was to develop a 
random terrain representation which could be replicated 
quickly but at the same time be representative of a particular 
"type" of terrain. The approach adopted for SIMTER generated 
a variable number of hill masses. The shape and height of 
the hills could be varied, but each hill was symmetric, since 
it was generated directly from the modified bivariate normal 
(MBVN) density function. 

SIMTER has two parts. The first part is the main 
program which creates terrain and, at the discretion of the 
user, plots both a three-dimensional drawing and a contour 
map. Used strictly for terrain generation, it can be used 
as a preprocessor for other combat models by producing grid 
points and elevations similar to those provided on computer 

^Needels , Parameterization of Terrain in Army Combat 
Analysis ♦ 



tapes by WES. The second part of SIMTER is a movement/line 
of sight routine which moves a target along specified routes 
across the generated terrain. 

C. A NEW REPRESENTATION METHODOLOGY 

The methodology developed in this thesis provides for 
the generation of hill masses utilizing MBVN density 
function to represent masses that are not symmetric. This 
development provides a significant increase in flexibility 
to represent hill masses more realistically. The mathematical 
development of this methodology is presented in the next 
chapter, after which actual results from executing the model 
are given. 
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III. THE EXTENDED MODEL 



A. UNSYMMETRIC TERRAIN 

The unsymmetric slope of terrain can be generated by 
a cutting method using the MBVN as follows: 



Z = MH-EXP \ - 



xa-p ) 
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maximum height 
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hill 


Z : 


height of hill 






y , y : 
x y 


spread of hill 
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ellipse factor 


of 


hill 



The unsymmetric slope of the hill can be represented by the 
vertical value of the MBVN subtracted from the vertical value 
of a linear function as shown in Figure 1 by the value HI. 




FIGURE 1. UNSYMMETRIC SLOPE GENERATION 
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Two basic questions must be addressed: 

1. Where do we start the cutting of the MBVN function? 

2. How do we compute the initial cutting position? 

If M is the direction of the cutting line (Line A) , then 
the cutting plane is through the location of the center of 
the hill as shown in Figure 2. 




FIGURE 2. CUTTING LINE DETERMINATION 

If the initial point of cutting is (XIN,YIN) and HI is 
height at the initial point, this point is on the cutting 
plane and also on the MBVN distribution function. Thus, if 
we consider the graph in two dimensions as shown in Figure 2, 
Line B is tangent to the circle at the point that the vertical 
value on the MBVN function is HI and the slope of Line B 
is -1/M. If the initial point of the cutting is (XIN,YIN) 
this point is on Line B, Line A and the circle at a height 
HI on the MBVN function. We can compute the initial point 
(XIN , YIN) as follows. 
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If the hill is cut from left to right and a = a , then 

x y 



XIN = 



u - a / ~ 2 ln(HI/MH) 
X X ' 1 + M 2 



( 8 ) 



YIN = y - a M 

y y 



-2 In (HI/MH) 



1 + m' 



(9) 



If the hill is cut from right to left and = a , then 



XIN = 



z2_ln(HI Affl) 
X X 1 1 + M 2 



( 10 ) 
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y + a M 

y y 



-2 In (HI/MH) 
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Thus, if the direction of cutting and HI are selected, then 
the starting point can be computed. Care must be taken in 
the selection of the value of HI. If HI is too small, then 
the hill will have an unrealistic shape. The value of HI 
affects the shape of the hill as shown in Figure 3. 




FIGURE 3. EFFECT OF HI ON THE HILL SHAPE 



15 



2. o < P < 1 Case 



HI = MH * EXP{ - 



2(l-p z ) 



■[(- 



XIN-y XIN-y YIN-y 

-) Z - 2 (— £) ( — £) 



X 



YIN-y 

+ (— — -Z) 



1 



(12) 



, XIN-y 9 

-2(1 - p ; In (HI/MH) = (— £) Z - 2 ( 



x 



XIN-y YIN-y 






x y 

YIN— y 

+ ( - ■ -■ V 



(13) 



From (13) the derivative is computed as follows: 
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FIGURE 4. CUTTING DIRECTION 

If the cutting direction goes through the longer 
axis as shown in Figure 4 , then 



cot 2a = 



_ 2 2 
o — a 

* — JL- 

-2po a 
x y 



( 22 ) 



The cutting direction, a, is determined by 



n 2 n 2 

, a - a 
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2 -2po 
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If the cutting direction goes through the shorter 
axis as shown in Figure 4, then 



6 = 90° - a 



(24) 
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Next, the height of intersecting hill masses must 
be determined as illustrated in Figure 5. 




XA 

FIGURE 5. INTERSECTING HILL MASSES 



HB defines height of hill "B" after cutting. 

HA defines height of hill "A" after cutting. 

Height of hill "A" is greater than that of hill "B" 
at the point XA but height of hill "B" is greater than that 
of hill "A" after cutting. If there is more than one hill 
located close to each other, those hills must be adjusted 
by comparing the height after cutting each one . The 
computation of the height of the terrain at any point is 
described in the following section. 

B. TERRAIN HEIGHT 

The height of the hill which is generated by the MBVN 
distribution corresponds to the height of a hill which is 
generated by the cutting method. The height of a hill 
at any point is computed as follows: 
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The height of the MBVN density function is given by 



= MH*EXP {- 
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The height of the linear function is 



(25) 



ZCUT = (X - XIN) tan 0 + (Y - YIN) tan <j> 



(26) 



where 6 = cutting angle of the x-axis 

4> = cutting angle of the y-axis. 

The height of the hill after cutting is 
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greatest height after cutting is- used. 



C. MOVEMENT LINE OF SIGHT 

The procedure of computing the line-of-sight between 
any two points on the terrain is described in this section. 
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If the position of the observer is (OX,OY) and the position 
of the target is (X,Y) , then the slope of the Observer- 
Target (0-T) Line is given by 



ML 



X - OX 
Y - OY 



(28) 



If the intersection point of a hill and the 0-T Line is 
(XI / YI) , then HI and the 0-T Line are perpendicular as 
shown in Figure 6 . 







FIGURE 6 . INTERSECTION OF 0-T LINE 
WITH HILL MASS 



Thus , 
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YI 



ML 2 • OX - ML • OY + y + ML 
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1 + ML 2 
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(29) 

(30) 
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If there is more than one hill on the 0-T Line, the 
dominant hill (after cutting) is used to compute the percent 
of the target covered relative to the observer- Target 
movement along the surface slope of the hill is computed 
as follows: 



ZH 



mh*exp { [(— £) 

2 (1-p J 0 x 




(X-XIN) tan 9 - (Y-YIN) tan <j> 
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If DEL is the direction of target movement, then 



DEL 



9ZH 

3x 



cos r 



+ 



3ZH 

ay 



sin 



r 



If the speed of the target is constant "V" 
then 



velocity = - 

I 



V 

1 + (DEL) 2 



If the speed of the target is dependent on 
the surface, then 

velocity = V + W (DEL) 

where W is the speed factor of slope. 



( 34 ) 



on the surface. 



( 35 ) 



the slope of 



( 36 ) 
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IV. TEST AND EVALUATION OF MODEL RESULTS 



A. DESIGN AND CUTTING ANGLE METHODOLOGY 

The methodology for parametrically describing terrain 
as a function of the number of hills, the center and spread 

of each hill, and the slope was developed in the previous 

chapters. In this chapter the effect of the number of hills, 
height of hills, and the slope on various intervisibility 
measures is investigated. Three random terrains were 
generated (with each terrain having six, twelve, and eighteen 
hills) resulting in nine basic terrain configurations. 

For each configuration three average hill heights (100, 

200, and 300 meters) were utilized. 

Finally, for each of the 27 resulting configurations, 
four slope cutting methods were used as follows: 

1. Down slope-cut from observer side to target side. 

2. Side slope-cut from right to left as viewed from 

the target. 

3. Up slope — cut from target side to observer side. 

4. Both slope-symmetric slope as viewed from any 
position. 

A total of 108 unique terrain configurations were 
generated for the investigation. For each of the config- 
urations, the mean and standard deviation describing the 
center and spread of each hill is specified. In addition, 
the cutting angle for determining slope must also be 
specified. 
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If the cutting angle is constant, then the hill can be 
described as follows: 

HILL = F(Peak, Position, Standard Deviation, Slope) 

In the constant cutting angle case, only the slope is 
variable, with peak, standard deviation, and position 
constant. In other words, if the slope varies in Figure 7, 
the standard deviation of the resulting distribution also 
varies . 




FIGURE 7. CUTTING ANGLE AND SLOPE 

The objective of the methodology is to maintain a constant 
standard deviation of the distribuition function . resulting 
from each cutting angle, so that only the slope of the hill 
varies. Otherwise the determination of how slope affects 
the intervisibility parameters would not be possible, since 
it would be confound with the standard deviation (spread) 
of the hills. 
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The problem, therefore, was to determine a method to 
use various cutting angles and maintain a constant standard 
deviation. After extensive investigations into various 
approaches to this problem, no method was found which would 
mathematically guarantee a constant standard deviation. 

The method shown in Figure 8 was adapted as the best 
approximation method available. 




In Figure 8, HT is the maximum height of the hill after 
cutting with Line A. Note that the cutting angle is now 
specified in reference to the center line of the original 
bivariate normal hill mass distribution. In order to 
generate a symmetric hill mass with approximately the same 
standard deviation (Line C) the cut is made maintaining 
a constant maximum height, HT, as in Figure 8. 

Several test runs made using this methodology indicated 
that the resulting standard deviations were very nearly 
equal for cutting angles. However, very little change in 
slope results for cutting angles greater than twenty degrees. 
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B. OBSERVER-TARGET POSITIONS 

For each terrain an observer position and initial target 
position were selected. Two routes from the target to the 
observer position were selected in accordance with accepted 
tactics of advance routes. Selected terrains indicating 
the observer-target position and selected routes are shown 
in Figures 12 through 17. in order to assure that the 
"number of hills" factor was properly analyzed, the same 
routes of approach on each generated terrain were used, 
independent of the number of hills. 



C . TARGET MOVEMENT 

Two cases of target movement were considered in the 
analysis . 

1. Constant Speed 

The target moves along the route at a constant speed 
of six meters per second, independent of hill slopes. The 
constant speed runs serve as the base case for measurement 
of the percent of time intervisible along the route. 

2 . Variable Speed 

The speed of the target is a function of the terrain 
slope being traversed. The speed of the vehicle is computed 
by the relationship as follows: 

velocity = 6.0 + 4.0 (slope) 

This slope is the directional derivative value at the 
instantaneous target position on the hill. 
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D. TARGET INTERVISIBILITY 



Five dependent variables describing target intervisibility 
are computed for each advance route as follows: 

1. Percent of time intervisible along the route, 

2. The number of intervisibility segments along the route, 

3. The average intervisibility segment length, 

4. The average intervisibility segment time, and 

5. The total length of intervisibility over the route. 
Intervisibility parameters are computed by utilizing the 
line of sight algorithm at fixed time intervals, TL, during 
the movement run. For this analysis, TL was set equal to 
one second. A sample output for one compute movement trace 
is given in Figure 9. 

E . EXPERIMENTAL DESIGN 

The objective of the analysis is to determine the effects 
of the number of hills (A) , the height of hills (B) and the 
slope of the hills (R) on target intervisibility. To 
accomplish this objective, an experiment was designed as 
described by the following design matrix. 

SLOPE 

# of Hills Height Down Side Up Both 

100 

6 200 

300 

100 

12 200 

300 

100 

18 200 

300 
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FIGURE 9. OUTPUT OF DATA 
Definition of Variables 



N 


Number of hills 


XMEL, YMEL 


Spread of hills 


HITE 


Maximum height of hill 


RUFX , RUFY , RUFP 


Different amount of the spread and height 
among hills 


INCX, INCY 


Grid increment for the terrain matrix 


XMAP , YMAP 


Size of the grid system in meters 


XGRID , YGRID 


Center position of hill 


HIGHT 


Maximum height of hill in modified bivariate 
normal distribution 


XSTD, YSTD 


Spread of hill (standard deviation) 


XSLOP , YSLOP 


Cutting angle 


DIRECTION 


Cutting direction 


XINIT, YINIT 


Initial position of cutting plane 


HEIGHT 


Hieght of hill at the center position on MBVN 
after cutting 


RX, RY 


Coordinates of approaching route of target 


TIME 


Movement time 


MAP COORDINATES Target position at given time 


ELEVATION 


Target elevation on the terrain 


SPEED 


Velocity on the terrain at the given time 


VISIBLE - 


1 means visible and 0 means invisible 



The different height of O-T line at the 
maximum height of surface of the terrain on 
the 0-T line, and the lower height between 
observer and target 



HLIN 



DZ 

KL 

DLIN 



Maximum height of surface of the terrain on the 
0-T line 

The hill which is maximum height between the 
observer and target 

Difference maximum height of surface between 
observer and target, and lower height between 
observer and target. 
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FIGURE 9 CONTINUED 
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The appropriate linear model is as follows : 



Y. . 



ijk - U + + R k + AB. j + AR. k + BR jk + ABR^ 



jk 



+ e. „ 
1 jk£ 



where 

Y ijkJl = i nterv i sibil ity dependent variable measure 

y = grand mean of the population 

A = number of hills, where 

1-6 hills 
i = 2-12 hills 

3-18 hills 

B = average height of hills, where 

1-100 meters 
j = 2-200 meters 

3- 300 meters 

R = direction of slope of hills, where 

1-down slope 
. _ 2-side slope 
" 3-up slope 

4- both (symmetric) 

Within each cell of the matrix, there are six replications 
consisting of two routes on each of the three basic terrain 
configurations . 

For the first treatment levels of slope defined above, 
a cutting angle of fifteen degrees was used, the symmetric 
case was generated with a cutting angle of zero degrees 
using the methodology previously described. Four analyses 
were conducted as representative of the large number of 
possible sitatuations . The measurements selected were 

1. Percent of time intervisible for constant speed 
and for variable speed, and 

2. Number of intervisibility segments for constant ^jpeed 
and for variable speed. 
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The results for each selected analysis are described in 
the following section. 



F. ANALYSIS OF RESULTS 

1 • Percentage of In ter visibility Time 
a. Constant Speed 



(x 



1/8 



Transgeneration Data) 

SUMMARY TABLE - ANOVA 



Source 


DF 


SS 


MS 


F Ratio 


# of hills 


2 


0.0437 


0.0218 


23.6417 


A 


height of hill 


2 


0.0117 


0.0058 


6.3209 


B 


slope of hill 


3 


0.0387 


0.0129 


13.962 


R 


AB interaction 


4 


0.0006 


0.0081 


0.1527 


AR interaction 


6 


0.0026 


0.0004 


0.4759 


BR interaction 


6 


0.0016 


0.00027 


0.2943 


ABR interaction 


12 


0.0012 


0.0001 


0.1114 


Error 


179 


0.1655 


0.0009 





Each of the three main effects above are significant 
at the 0.01 level, but no interaction terms are significant. 

The following table gives the mean and standard deviation 
for each treatment level of slope. 
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Slope 



Mean 



Std. Dev 



down 


0.6130 


0.1539 


side 


0.5214 


0.1701 


up 


0.6926 


0.1629 


both 


0.5573 


0.1595 



Note that the up slope has the greatest percentage of inter- 
visibility time; down slope, both, and side slope following 
in that order. 

The following table gives the mean and standard 
deviation for the treatment levels of number of hills and 
height of hills for the constant speed case. 



# of Hills 


Mean 


Standard Deviation 


6 


0.6976 


0.1620 


12 


0.5462 


0.1620 


18 


0.5272 


0.1490 


Height of Hill 


Mean 


Standard Deviation 


100 


0.6475 


0.1775 


200 


0.5844 


0.1684 


300 


0.5556 


0.1628 
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b. Speed which Depends on Slope 



, 3/8 _ 

(x Transgeneration Data) 



SUMMARY TABLE - AN OVA 



Source 


DF 


SS 


MS 


F Ratio 


# of hills 


2 


0.5152 


0.2576 


34.286 


A 


Height of hill 


2 


0.2041 


0.1020 


13.582 


B 


Slope of hill 


3 


0.3830 


0.1277 


16.991 


R 


AB interaction 


4 


0.0033 


0.0008 


0.1090 


AR interaction 


6 


0.0381 


0.0064 


0.8450 


BR interaction 


6 


0.0165 


0.0028 


0.3670 


ABR interaction 


12 


0.0091 


0.0008 


0.1000 


Error 


179 


1.3449 


0.0075 





Each of the three main effects above are significant at the 
0.01 level, but no interaction terms are significant. 

c. Constant and Dependent Speed!. 

The ANOVA for comparing constant and dependent 
speed is given in Figure 10. The treatment level data is 
as follows: 

Mean Standard Deviation 

Constant Speed 0.569 0.1732 

Dependent Speed 0.5393 0.1921 
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2 . Number of Intervisibility Segments 



a. Constant Speed 

SUMMARY TABLE - ANOVA 

1/2 

(X Transgeneration Data) 



Source 


DF 


SS 


MS 


F-Ratio 


# of hill 


2 


8.3415 


4.1708 


82.281 


A 


Height of hill 


2 


0.2402 


0.1201 


2.369 


B 


Slope of hill 


3 


2.7278 


0.9093 


17.938 


R 


AB interaction 


4 


0.0659 


0.0165 


0.3255 


AR interaction 


6 


0.2686 


0.0448 


0.8838 


ABR interaction 


12 


0.1389 


0.0116 


0.2288 


Error 


179 


9.0475 


0.5056 




From the ANOVA, the 


number 


of hills and 


slope are 


significant 


at the 0.01 level, 1 


but the 


height of hills is not 


significant 
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b. Dependent Speed 



• SUMMARY TABLE - AN OVA 

1/2 

(X Transgeneration Data) 



Source 


DF 


SS 


MS 


F-Ratio 


# of hill 


2 


8.2224 


4.1112 


79.337 


A 


Height of hill 


2 


0.2561 


0.1281 


2.4272 


B 


Slope of hill 


3 


2.6384 


0.8795 


16.972 


R 


AB interaction 


4 


0.06446 


0.01611 


0.311 


AR interaction 


6 


0.2905 


0.0484 


0.934 


ABR interaction 


12 


0.1289 


0.0107 


0.2065 


Error 


179 


9.2756 


0.0518 


0.2065 


The results of the dependent 


speed case 


are analogous to 


those of the constant 


speed 


case . That 


is, the 


number of 



hills and slope of hills are significant, but the height of 
hills is not significant . 

The mean and standard deviation of the treatment levels 
are given below. 
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SLOPE 



MEAN 



STD DEV 





CONST 


DEP 


CONST 


DEP 


DOWN 


4.2075 


4.2075 


1.04444 


1.0444 


SIDE 


4.1852 


4.1852 


1.1339 


1.1339 


UP 


3.3148 


3.3333 


1.1947 


1.2131 


BOTH 


4.4259 


4.4259 


1.3818 


1.3818 


# OF HILL 


MEAN 

CONST 


DEP 


STD 

CONST 


DEV 

DEP 


6 


4.930 


3.038 


1.1871 


0.7166 


12 


4.139 


4.139 


1.039 


1.0388 


18 


3.042 


4.931 


0.7207 


1.1788 


HEIGHT OF 
HILL 


MEAN 

CONST 


DEP 


STD 

CONST 


DEV 

DEP 


100 


3.875 


3.875 


1.393 


1.393 


200 


4.069 


4.069 


1.191 


1.191 


300 


4.155 


4.169 


1.191 


1.195 
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V. CONCLUSION AND RECOMMENDATION 



Terrain is one of the key factors in the analysis of 
land combat. Thus, the representation of terrain in a land 
combat model is a critical factor in the model design. High 
resolution combat models typically employ actual terrain 
data to describe the terrain environment. This approach, 
although realistic for the specific terrain being represented, 
is very costly and time consuming to implement. 

The goal of the research by Needels and this thesis is 
to provide a parametric representation of terrain usable 
either as input to a high resolution combat model or as a 
"stand-alone" model to investigate intervisibility segment 
length parameters. In particular, the current model 
describes the terrain as a function of the number of hill 
masses, the height of the hills, the position and spread 
of the hills, and the slope. 

A MBVN distribution function and a linear function 
cutting the MBVN function are employed to achieve the 
terrain representation. The terrains given in Figures 12 
through 17 illustrate the flexibility of the model. In 
fact, with appropriate selection of the input parameters 
(given in Figure 9) any desired terrain may be generated. 

The results of the analysis described in Chapter IV 
indicate that slope is highly significant in the resulting 
values of intervisibility parameters . Some of the trends 
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noted from the results given in Chapter IV are described 
below. 

1. An up slope (cut from target side to observer side) 
resulted in the largest value of percent intervisibility 
time . 

2. Increases in the number and height of hills resulted 
in decreases in the percent intervisibility time. 

3 . Speed dependent on terrain slope resulted in a 
slight reduction in percent intervisibility time from the 
constant speed case. 

The analysis of the hypothetical data described in 
Chapter IV serves to illustrate the potential utilization 
of the model. The trends indicated from the analysis are 
appropriate for the selected terrains and conditions, but 
should not be construed as generalized results. 

An alternate method (tilted method) for generating 
unsymmetric terrain outlined below is suggested as an area 
for future enhancement of the model. Given a MBVN density 
function centered at Y in Figure 11, define 0 as the tilted 
angle, assuming p = 0 . The height of the density function 
is Z2 at Y2 and is given by 




(38) 



and 



Z2 



tan - 0] 



(39) 



Y2-Y 



4-2 




FIGURE 11 TILTED METHOD 



From 38 and 39, 

MH*EXP {- | [(-^V + (-^n) = (y-y 2 ) tan ( J - e) (40) 

For a specified value of x and y, the value of y 2 can be 
computed using the golden section method. The tilted height, 
ZH, can then be computed as follows: 



ZH 



sin 9 



(41) 



Additional investigation of this method is required before 
implementation in the model. 

In conclusion, a model has been developed and exercised 
which provides a potentially valuable tool for efficient 
analysis of various terrain effects in the land combat 
environment. The intervisibility parameters output from 
the model can serve as critical inputs to land combat models 
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for evaluating such things as the contribution of mobility, 
agibility, and firepower to battlefield survivability. 
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FIGURE 12-2. Terrain: 1, # of Hills; 18, Average Height 
100, Slope; Dqwh 
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Height; 300, Slope; Up 
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Terrain: 2, # of Hills? 18, Average Height; 
300, Slope; -Up 
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FIGURE 13-2. 
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FIGURE 14-2. Terrain: 2, # of Hills; 12, Average Height; 
100, Slope; Up 
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FIGURE 16-2. Terrain: 3 # of Hill?; 6_, Average Height; 

300 , Slope; - Side.^ 
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